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Quantum entanglement under an extensive bipartition can reveal the critical boundary theory of a topological 
phase in the parameter space. In this study we demonstrate that the infinite-randomness fixed point for spin-1/2 
degrees of freedom can emerge from an extensive random bipartition of the spin-1 Affleck-Kennedy-Lieb-Tasaki 
chain. The nested entanglement entropy of the ground state of the reduced density matrix exhibits a logarithmic 
scaling with an effective central charge c = 0.72 ± 0.02 « In 2. We further discuss, in the language of bulk 
quantum entanglement, how to understand all phase boundaries and the surrounding Griffiths phases for the 
antiferromagnetic Heisenberg spin-1 chain with quenched disorder and dimerization. 
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I. INTRODUCTION 

Entanglement spectrum under an extensive bipartition of 
a topological ground state has recently emerged as a novel 
approach to study the quantum phase transition between the 
topological phase and its trivial counterpartjivi The so-called 
bulk entanglement spectrum (BES) reveals the boundary the¬ 
ory in the corresponding parameter space of a model system, 
rather than the edge states along the fictitious boundary in real 
spacei2ii£ These studies suggest that a triangular correspon¬ 
dence among the bulk theory, the edge theory, and the critical 
theory may exist generically for a topological phase.-^iii 

An instructive example is the spin-1 Haldane gapped 
phase,whose fixed-point properties are captured by the 
valence-bond solid picture of the Affleck-Kennedy-Lieb- 
Tasaki (AKLT) model Hamiltonian.'^ The corresponding ex¬ 
act ground state wave function can be expressed as a matrix 
product state (MRS), indicating that the fundamental entities 
of the spin-1 chain is actually the fractionalized spinons carry¬ 
ing spin 1/2 that form singlets across adjacent sites. The spin- 
1/2 object can be observed as the edge mode at the end of the 
chain, or at the end of a partition in the entanglement study. 
The degeneracy due to the fractionalized spin in the entan¬ 
glement spectrum of the antiferromagnetic Heisenberg spin-1 
chain is key to understand the Haldane phase as a symme¬ 
try protected topological phasejii Under a uniform extensive 
bipartition with disjoint segments, the segment end spins co¬ 
alesce into an emergent critical spin-1/2 chain with a central 
charge c = Ir which describes the collapse of the Haldane 
phase under imposed dimerization^ (marked by the stars in 

Fig IB- 

Meanwhile, spin chains can be interesting with a random 
probability distribution of the bond couplings between neigh¬ 
boring spins, which may be broaden without limit as the sys¬ 
tem is coarse-grained.*^ Such a system is governed by the 
infinite-randomness fixed point (IREP), whose ground state 
is characterized by a random pattern of spin singlets formed 



FIG. 1. (Color online) Phase diagram of the antiferromagnetic 
Heisenberg spin-1 chain in the presence of quenched bond disorder 
and dimerization.^’ The regimes of interest for the extensive biparti¬ 
tions in Fig.l^are labeled accordingly. The red stars mark the critical 
spin-1/2 chain that can be accessed by the extensive uniform biparti¬ 
tion of the AKLT state.^ The thick black boundaries are the random 
singlet states that can be accessed by extensive random bipartitions 
of the AKLT state (see text). 

over large spatial separationsJ^^ For the antiferromagnetic 
Heisenberg spin-1 chain the Haldane phase is stable against 
weak disorderi^iiS and weak dimerizationi^S but the emer¬ 
gent critical spin-1/2 chain is not»i^iii As illustrated in Fig.[T] 
the Haldane phase and two dimer phases in the presence of 
quenched disorder are separated by the spin-1/2 random sin¬ 
glet (RSiy 2 ) boundaries, which merges into a single spin-1 
random singlet (RSi) boundary between the dimer phases. 
In addition, unusual Griffiths effects, characterized by two 
continuously varying dynamical exponents, appear near the 

boundaries 

In this paper we will show that disorder physics can also 
be revealed in the bulk entanglement study under an extensive 
random bipartition of the spin-1 chain. For uncorrelated seg- 
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ment length, the nearest-neighbor couplings in the bulk entan¬ 
glement Hamiltonian exhibits a power-law distribution, which 
is precisely the fixed-point solution for random antiferromag¬ 
netic spin-1/2 chains under strong-disorder renormalization- 
group (SDRG) transformation.'^ Therefore, the ground state 
of the entanglement Hamiltonian realizes the RS 1/2 state, for 
which we provide further evidence by fitting the ensemble- 
averaged nested entanglement entropy to a logarithmic scal¬ 
ing with an effective central charge c = 0.72 ± 0.02 In 2. 
We discuss how to vary the random bipartition to explore 
the phase diagram of the antiferromagnetic Heisenberg spin-1 
chain with quenched disorder and dimerization. 

II. THE AFFLECK-KENNEDY-LIEB-TASAKISPIN-I 
CHAIN 

The spin-1 AKLT parent Hamiltonian on a periodic chain 
of length L is defined 
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whose exact ground state can be expressed as an MRS 
I'Raklt) = ^ Tr • • • aI^^I) |siS 2 • • • sl) , 


Si • Si_|_i + — (Si • Si_|_i) 


( 1 ) 


for J > 0, where the local physical spin Si = —1,0, -1-1 and 
A[®1 are local 2x2 matrices given, e.g., in Ref.S In the ther¬ 
modynamic limit, the spin-spin correlation function decays 
exponentially with a correlation length ^ = l/ln3 sa 0.91, 
and any two spins that are separated by an even number of the 
lattice sites are antiferromagnetically correlated. 

To study the quantum entanglement of the AKLT state, we 
can divide the chain into A and B partitions, and define the 
entanglement Hamiltonian He through the reduced density 
matrix 


PA = Trs (|T'aklt)('1'aklt|) = e (2) 

In a common practice, as illustrated in Fig. |^, A is an open 
spin segment of length I > 2. The MRS representation dic¬ 
tates that PA contains four eigenvalues: a singlet and a triplet, 
as required by the SU(2) symmetry. Explicitly, the singlet and 
triplet eigenvalues ar e^^'^^ 



The corresponding entanglement Hamiltonian can thus be de¬ 
termined to be He = ■ te, with 


J{1) = In 


Li-HH 
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FIG. 2. (Color online) Various bipartitions of the AKLT state and 
the resulting effective spin models for the entanglement Hamilto¬ 
nian. We denote A the collection of segments with orange shadow 
and B that of the segments with green shadow, (a) The entanglement 
spectrum (illustrated in the box) of a finite segment with an even 
number of sites is equivalent to the spectrum of two antiferromag¬ 
netically coupled spin-l/2s. (b) An extensive uniform bipartition of 
the AKLT chain (with even-length segments) generates a pure anti- 
ferromagnetic spin-1/2 chain, whose ground state can be described 
by a c = 1 CFT. (c) An extensive random bipartition with even- 
length-only segments generates an effective random antiferromag¬ 
netic spin-1/2 chain, whose ground state is the RS 1/2 state, and (d) 
the corresponding pA in the MRS representation, (e) An extensive 
random bipartition with isolated spins separated by segments of even 
length generates an effective random antiferromagnetic spin-1 chain, 
whose ground state is the RSi state. 


where we have taken the large I limit, and introduced Jq = 4 
as the energy unit, ^ = 1/ In 3 is the correlation length, 
and te are the two fractionalized spin-l/2s at segment ends. 
Hence, they are coupled antiferromagnetically for even I, and 
ferromagnetically for odd 1. If we swap A and B, we obtain 
the identical entanglement Hamiltonian, which means that the 
coupling between the two end spins is independent of their 
being connected by a segment in A or by that in B. 

On the other hand, the bulk entanglement Hamiltonian He 
stems from the extensive bipartitions that divides the spin 
chain into A and B sets of alternating segments (see Fig. |2j) 
for an example). The word “bulk” emphasizes that the bound¬ 
aries between A and B spread out the whole system, and the 
couplings between the boundary spin-1/2s are relevant. Due 
to the gap of AKLT model, the coupling strength between 
spin-l/2s decreases exponentially with their distance, hence 
only nearest neighboring interactions dominate in large I limit, 
and the bulk entanglement Hamiltonian is described by the 
spin-1/2 Heisenberg model, with coupling strength given in 
Eq. (H.s Hence, the He with even segment length is critical 
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and governed by an SU (2)i Wess-Zumino-Witten conformal 
field theory (CFT) with central charge c = 1. This theory is, 
however, unstable against arbitrarily small disorder. 


III. RANDOM BIPARTITIONING OF THE AKLT CHAIN 

A naive thinking to study the effect of quenched disorder 
is to randomize the bond couplings in the parent Hamiltonian, 
such that we modify the ground state wave function, hence 
the reduced density matrix. The AKLT case is, however, an 
exception. The entanglement spectrum is immune to weak 
bond randomness [i.e., if we replace uniform J with random 
Ji in Eq. ([T]i] as the exact ground state of the random-bond 
AKLT model is identical to that of the pure model.This is a 
vivid example that the perturbation can influence thermal ex¬ 
citations but not the excitations in the entanglement spectrum. 

This motivates us to enforce disorder by introducing an ex¬ 
tensive random bipartition for the AKLT state, as illustrated in 
Fig.|2};. For the initial simplicity, we assume that the number 
of lattice sites in each segment is even. After tracing out the 
degrees of freedom in every other segments, the bulk entan¬ 
glement Hamiltonian now reduces to a spin-l/2s Heisenberg 
model with random couplings, that is 

He ■ Ti+i, (5) 

i 

where Ji = J {k) with li being the length of the ith seg¬ 
ment, and the form of J (1) is given in Eq. (|4|i. On physi¬ 
cal ground we may then expect that the low-energy physics is 
governed by the IRFP. Normally, this fixed point is revealed 
by a real-space decimation process, commonly referred to as 
the SDRG, developed by Dasgupta and Mai^ and by Fisher. 

In general, the SDRG process in the initial stage depends 
strongly on the probability distribution of the random cou¬ 
plings, thus often drives the energy scale of interest to be ex¬ 
ponentially small than the strongest bond in the bare Hamil¬ 
tonian.^'* Surprisingly, this is not the case here because of the 
uncorrelated locations of the segment ends (which we impose) 
and the finite correlation length in the topological phase. 

For the distribution of segment length (hence the coupling 
strength), the most natural recipe is to specify a fixed aver¬ 
age length of the segments and to assume that the segment 
ends are located independently. If we suppose that the av¬ 
erage length of the segment is large, one can show that the 
probability distribution of the segment length I (assumed to 
be continuous for simplicity) satisfies the following differen¬ 
tial equation 

ip{t) pm 

where the constant I is the average length of the segments. We 
can further assume that the segments have a minimum length 
of Iq, such that 

P{1) ~ (7) 

The discreteness of the segment length is not important when 
the average segment length is sufficiently long. As often 


emerged in the minimum mathematical model of the life ex¬ 
pectancy problem, the exponential form of the probability dis¬ 
tribution simply means that the occurrence of the next seg¬ 
ment end is independent of the location of the previous one. 

As we discussed above, the energy scale that couples two 
adjacent segment end spin-1/2s depends on the length of the 
segment. When the segment length I (measured by the number 
of sites) is sufficiently long the effective coupling asymptot¬ 
ically approaches Eq. (0. For even I the effective couplings 
are all antiferromagnetic. We point out, though, including fer¬ 
romagnetic bonds can lead to a different SDRG fixed point 
with large spins formed in a random-walk fashion^ 

The length scales in the two previous exponential laws 
should in general be different. The former is the average seg¬ 
ment length enforced externally, while the latter is the inter¬ 
nal correlation length associated with the topological phase. 
Together, the probability distribution of the effective nearest- 
neighbor couplings reads 


P(J) 



( 8 ) 


where F = I/^ and H Once again we stress that, 

in Eq. (|5]l, the distribution of couplings Ji is random and is 
given by Eq. ®. 

One immediately recognizes that Eq. ® bears a marked 
resemblance to the fixed point distribution that describes the 
random singlet phase in the disordered S — 1/2 Heisenberg 
antiferromagnet, derived by Fisher-i^^ by using a real space 
based SDRG scheme. 

In Fisher’ s'll*^ SDRG scheme the fixed point distribution 
that characterizes the random-singlet phase is obtained as the 
solution of the master equation that describes the flow of the 
distribution function of coupling constant; 

5P(/3) _ 5P(/3) 
dr ~ d/3 

pOO pOO 

+ P(0) / d/3i / dp2P{/3i)P{/32)50,+p,-^, 

Jo Jo 

(9) 


where we follow Fisher to introduce the dimensionless scal¬ 
ing variable P = In(r2/J). The differential equation states 
that during the bond decimation the flow of the bond distri¬ 
bution has two contributions: (i) a shift in /3 due to the re¬ 
duction of the UV cutoff H, and (ii) the replacement of deci¬ 
mated bonds by the effective couplings generated through the 
second-order perturbation theory as the logarithmic RG flow 
parameter F changes. As alluded to earlier the Eq. (|9]l admits 
a fixed point solution of the form given in Eq. ®. However, at 
this point it is important to emphasize that this striking result 
elucidated in Eq. ®, is not a result of any SDRG calcula¬ 
tion that we have performed but due to the random bipartition 
scheme that we have adopted wherein the segment lengths are 
controlled by the distribution given by Eq. (|7]). We will com¬ 
ment in detail about the apparent similarity of Eq. 0, and the 
fixed point distribution of the random singlet phase derived 
by Fisher later on in the manuscript. At this juncture, it suf¬ 
fices to re-emphasize that the extensive bipartition results in 
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an entanglement Hamiltonian that corresponds to an S' = 1/2 
disordered Heisenberg model. 

The ground state entanglement entropy of the random 
Hesienberg model of Eq. (|5ll was calculated by Refael and 
Moora^^ by laying recourse to Fisher’s formulatioi>ii of the 
SDRG.*® They realized that the entanglement of a segment 
of size £ with rest of the infinite random Heisenberg chain 
is just the product of the number of singlet bonds spanning 
across the boundary times the entanglement entropy per sin¬ 
glet (which is In 2). To accurately determine the proportion¬ 
ality constant the history of the singlet bond formation across 
the boundary separating the two subsystems should be kept 
in detail. In Ref. |2^ the SDRG methodology was very suc¬ 
cessfully adopted to do the same, and it was shown that akin 
to critical systems, the entanglement entropy in the random- 
singlet phase of the S' = 1/2 Heisenberg magnet exhibits a 
logarithmic scaling with respect to the sub-system size £. In 
other words. 


S{£) = —— ln£ + constant. (10) 

o 

Here, the coefficient In 2 that controls the logarithmic behav¬ 
ior is interpreted to be an effective central charge c. 

Now, in Sec. lYl we will apply the random bipartition 
scheme and numerically calculate the entanglement entropy 
of the ground state of the resultant entanglement Hamiltonian. 
In particular, we compute an effective central charge and com¬ 
pare it to that of the disordered Heisenberg chain c = In 2 [see 
Eq. (ITOl) above]. 


IV. NUMERICAL PROOF 

To demonstrate that the IRFP is indeed accessible even in 
small systems, we study the nested entanglement entropy, a 
straightforward generalization of the uniform case,^ and com¬ 
pare it to the entanglement entropy of random S = 1/2 
Heisenberg chains, which is controlled by an effective cen¬ 
tral charge c = ln2i^ For this purpose, we impose peri¬ 
odic boundary condition and compute pA, whose matrix prod¬ 
uct representation is illustrated in Fig. |2}l, then determine its 
ground state. We further divided the subsystem A with n seg¬ 
ments of a total La spin-I s into left partition (p segments) and 
right partition (n — p segments) 

^ ^ ^ ^ ^ ^ 

' = E?=1 li 

with I and La — I spins, respectively. Under this nested bi¬ 
partition of A, We calculate the nested entanglement entropy 
s{l, La) for the ground state of pA- We average over random 
realizations of the segment lengths in A and B, among samples 
with the same number of segments and the same pattern of the 
nested bipartition, and plot the averaged s{l, La) as a function 
of sample averaged Inp((, = ln[(L^/7r) sin(7rZ/L^)] in 

Fig. [3 for samples with 4-12 segments of random length pre¬ 
scribed by the discrete version of Eq. (|7]) with Iq = 2 and 
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FIG. 3. (Color online) The ensemble-averaged nested entangle¬ 
ment entropy s{l, La), as a function of the averaged \ng{l, La) = 
ln[(LA/7r) sin{nl/LA)] for the extensive random bipartition, where 
La is the size of A, while I and La — I are the sizes of the nested 
subsystems. The lengths of the segments are chosen according to 
the discrete version of Eq. 0 with lo = 2 and / = 5. We also 
plot the uniform case^ for comparison. The slopes of the linear fits 
are 0.242 ± 0.007 and 0.338 ± 0.004, corresponding to (effective) 
central charges c « In 2 and c « 1, respectively. 


1 = 5^^. Depending on segment number, we choose 2,000- 
10,000 random realizations. For comparison we also show 
data for the uniform bipartition with two spins in each seg¬ 
ment.^ The data in the random case can be fitted by 


s{1,La) = -In 


La . 

— sin 

TT 



( 12 ) 


where the effective central charge c = 0.72 ± 0.02 is in excel¬ 
lent agreement with the expected value of In 2 0.693, as op¬ 

pose to c = 1 in the uniform caset» The result is remarkable as 
the largest system contains a mere 12 segments with average 
segment length 1 = 5 and without involving SDRG. The nu¬ 
merical result also demonstrates that the exponentially small 
longer-range couplings among the segment end spin-1/2s are 
indeed irrelevant. 

The instant arrival at the fixed-point solution of the bulk en¬ 
tanglement Hamiltonian is very reminiscent to the case of the 
extensive uniform bipartition.^ For the concrete example in 
Fig.|2}r, one can imagine that a spin-1/2 singlet bond is formed 
within each segment, leaving an extra pair of spin-l/2s at the 
two ends to be coupled with the end spins in other segments, 
as the result of the bipartition. This maps the partition A (a 
spin-1 chain with bipartition-induced critical dimerization) to 
a spin-1/2 chain with an effective dimerization 5 = 0, since 
the lengths of A segments and B segments are equal. There¬ 
fore, we arrive instantly at the critical point for the pure spin-1 
chain with imposed dimerization (i.e., not due to spontaneous 
symmetry breaking). 

Now, in the case of the random bipartition the almost in¬ 
stantaneous arrival at the fixed point solution is also predicated 
on one another important condition: Namely, our choice of the 
probability distribution for the segment length [see Eq. (|7])]. It 
is this choice of probability distribution that ensures that our 
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FIG. 4. (Color online) The ensemble-averaged nested entangle¬ 
ment entropy s{l, La), as a function of the averaged lng{l, La) = 
ln[(LA/7r) sin(7r(/LA)] for binary partition, where La is the size of 
A, while I and La — I are the sizes of the nested subsystems. The 
length of each segment is chosen to be Li =4 with probability p 
or 1/2 = 6 with probability 1 — p. The effective central charge c 
can be extracted from the slope to be 0.75 ± 0.03 for p = 0.1 and 
0.81 ± 0.03 for p = 0.5. The values of c are in between c = 1 
(pure case) and c = In 2 (IRFP) as indicated by the two dashed lines 
(guide to eye only). 


results immediately converge to the fixed point distribution 
P{J) for the nearest-neighbor coupling elucidated in Eq. ®. 
Consequently, the random bipartition effectively sets up the 
bond distribution that looks akin to the fixed point distribution 
that characterizes the random singlet phase of the S — 1/2 
Heisenberg antiferromagnet. This provides a reason of why 
the numerically computed effective central charge c is almost 
bereft of finite size effects and locks on immediately to the 
universal value of In 2. Of course, this brings out immediately 
the question of how quickly one can arrive at the random sin¬ 
glet state fixed point if we were to start with a distribution 
that is non-Poissonian. We defer this discussion to the subse¬ 
quent paragraph. We end this paragraph with a few comments 
on the Poissonian distribution of segment length, Eq. Q. In 
our analysis we have ignored the fact that the segment length 
is discrete. This can be justified if we imagine approaching 
the fixed point via an SDRG: Under this SDRG, the (~ 1/ J) 
grows without limit, which means that the system flow to infi¬ 
nite randomness. In our setup this points our primary interest 
to the regime where I is large, hence our assumption that the 
discreteness of the segment length is not important and our 
simplification of the length dependence of the effective cou¬ 
plings J (1) between the segment end spins are justified. In the 
same limit the choice of Iq, which sets the UV cutoff, is also 
not important. 

To understand how the distance to the random singlet state 
depends on the segment length distribution, we consider a 
binary partition. The length of each segment is chosen to 
be Li with probability p or L 2 with probability 1 — p (with 
p < 0.5). Apparently, the binary distribution is not desirable 
as the initial distribution for an SDRG, as the probability that 
two neighboring lengths are equal can be high. Hence, the 
perturbative decimation is not valid. While it seems coun- 


p 

0.1 

0.2 

0.3 

0.4 

0.5 

c 

0.75 

0.75 

0.76 

0.77 

0.81 

A, 

0.03 

0.04 

0.03 

0.04 

0.03 


TABLE I. The effective central charge c with uncertainty Ac for 
binary partition. The length of each segment is chosen to be Li = 4 
with probability p or L 2 = 6 with probability 1 — p. 

terintuitive, one can argue that the binary distribution with a 
smaller p is closer to the IREP distribution than that with a 
larger p. The argument goes as follows; With a small p, we 
expect (with high probability) to have consecutive segments 
of length L 2 followed by a single segment of length Li, then 
another consecutive segments of length L 2 followed by a sin¬ 
gle segment of length Li, and so on. In the limit of small p, 
the number of consecutive length-L 2 segments N has a broad 
distribution (Poisson, indeed). Within a consecutive length- 
L 2 segments (say, of total length NL 2 ), we have effectively 
A^ -F 1 antiferromagnetically coupled segment-end spin-l/2s, 
with low-energy excitations being magnons with energy scale 
/N). Consequently, a broader distribution of en¬ 
ergy scales exists naturally in the entanglement Hamiltonian 
for smaller p. The argument works well for Li ^ L 2 , when 
consecutive segments of L 2 are coupled weakly. Neverthe¬ 
less, even for Li < L 2 , the conclusion is still valid, as the 
end spin-l/2s of the Li segments form singlets and generate 
much weaker couplings between their neighboring segment- 
end spins. In Pig. |4l we plot the ensemble-averaged nested 
entanglement entropy s{l, La), as a function of the averaged 
hig{l,LA) = \'n[{LA/T^) sln{'Kl/L a)] for the binary partition 
with p = 0.1 and 0.5. We choose Li = 4 < L 2 = 6 and find 
that the effective central charge c = 0.75 ± 0.03 forp = 0.1 
and c = 0.81 ± 0.03 for p = 0.5. We list the effective cen¬ 
tral charge c for five different p in Table |I] As expected, the 
smaller p is, the closer c is toward In 2 for the IREP value. 
On the other hand, the larger p is, the closer c is toward 1 for 
the disorder-free value. We note that, as we have only up to 
12 segments, finite-size artifacts can become significant for 
p<0.1soc = ln2 cannot be obtained in our calculation for 
binary partition. 

V. DISCUSSIONS AND CONCLUSIONS 

The entanglement-revealed criticality in both the pure case 
and the random case establishes that exploring the appropri¬ 
ate bulk entanglement Hamiltonian of the ground state wave 
function representing a topological phase can efficiently dis¬ 
till the quantum information belonging to the corresponding 
critical point separating the topological phase and its adjacent 
trivial phase. This is also observed for the integer quantum 
Hall transitions that can be classified by a Z index.* 

The distribution of the segment lengths can introduce addi¬ 
tional relevant perturbations for the emergent degrees of free¬ 
dom. One example is that the difference between the average 
segment lengths in A and B can introduce an effective dimer¬ 
ization 5 for the RS 1/2 state. Unlike the pure case, 5 can be 
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continuously tuned in the random case. The resulting bulk 
entanglement Hamiltonian effectively describes random spin- 
1/2 chains with weak but continuously varying bond dimer¬ 
ization. This can lead to critical behavior but with finite spin 
correlation length controled by 5, the characteristics of a Grif¬ 
fiths phase.— On the other hand, the random antiferromag¬ 
netic S' = 1 chains with enforced dimerization also contain 
spin-1 degrees of freedom, which, together with the spin-1/2 
degrees of freedom, result in Griffiths phases with two in¬ 
dependent dynamical exponents. The entanglement analogy, 
thus, requires that partition A (or B) should contain segments 
with odd length. We leave the details aside but point out that 
one can compare the extensive random bipartition case to the 
domain-wall description of Damle and Huse^^ hence the rel¬ 
evant physics will follow. 

It is interesting to point out that a spin-1 random singlet 
phase without dimerization^^ can also be accessed, e.g., by a 
bipartition (see Fig.|2^) that contains single sites in A (creating 
effective spin-1 s) and segments with random but even length 
in B (providing random antiferromagnetic couplings between 
effective spin-Is). Indeed, it can be shown the entanglement 
Hamiltonian is dominated by the antiferromagnetic Heisen¬ 
berg coupling between nearest neighboring spin-Is. This, 
however, does not fit the scenario that we discussed above, 
because the RSi state, as illustrated in Fig.[T] is a critical line 
that is separated from the Haldane phase (including the Grif¬ 


fiths region) by a multicritical point. The multicritical point 
has an emergent permutation symmetry corresponding to the 
interchange of the Haldane phase and the dimer phases that 
meet at the pointi^^ We leave it to future investigation whether 
it is possible to incorporate the permutation symmetry and to 
access the multicritical point. 

In summary, we have discussed how to understand the crit¬ 
ical phase boundaries of the antiferromagnetic Heisenberg 
spin-1 chain in the presence of quenched disorder and dimer¬ 
ization in terms of quantum entanglement by applying exten¬ 
sive bipartitions to the AKLT state that represents the topolog¬ 
ical Haldane phase. In particular, the calculation of the effec¬ 
tive central charge provides a striking example how efficient 
quantum entanglement can access the critical information of 
a topological quantum phase transition. 
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